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ABSTRACT 



^ \ Color-color diagrams for the clump and interclump emission in 10 dump- 

ed '. cluster galaxies of the Ultra Deep Field are made from B,V,i, and z images 



and compared with models to determine redshifts, star formation histories, and 



O ! galaxy masses. These galaxies are members of a class dominated by 5 to 10 giant 

^ \ clumps, and having no exponential disk or bulge. The redshifts are found to be 

in the range from 1.6 to 3. The clump emission is typically 40% of the total 
galaxy emission and the luminous clump mass is 19% of the total galaxy mass. 
O \ The clump colors suggest declining star formation over the last ~ 0.3 Gy, while 

\ the interclump emission is redder than the clumps, corresponding to a greater 

^ ' age. The clump luminous masses are typically 6 x 10® Mq and their diameters 

average 1.8 kpc, making their average density ~ 0.2 M© pc~^. Including the 
^ \ interclump populations, assumed to begin forming at z = 6, the total galaxy 

\ luminous masses average 6.5 x 10^° Mq and their diameters average 19 kpc to 

the 2 — 0" noise level. The expected galaxy rotation speeds average ~ 150 km 
s~^ if they are uniformly rotating disks. The ages of the clumps are longer than 
their internal dynamical times by a factor of ~ 8, so they are stable star clusters, 
but the clump densities are only ~ 10 times the limiting tidal densities, so they 
could be deformed by tidal forces. This is consistent with the observation that 
some clumps have tails. The clumps could form by gravitational instabilities in 
accreting disk gas and then disperse on a ~ 1 Gy time scale, building up the 
interclump disk emission, or they could be captured as gas-rich dwarf galaxies, 
flaring up with star formation at first and then dispersing. Support for this 
second possibility comes from the high abundance of nearly identical clumps in 
the UDF field, smaller than 6 pixels, whose distributions on color-magnitude and 
color-color plots are the same as the galaxy clumps studied here. The distribution 
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of axial ratios for the combined population of chain and clump-cluster galaxies 
in the UDF is compared with models and shown to be consistent with a thick 
disk geometry. If these galaxies evolve into today's disk galaxies, then we are 
observing a stage where accretion and star formation arc extremely clumpy and 
the resulting high velocity dispersions form thick-disks. Several clump-clusters 
have disk densities that are much larger than in local disks, however, suggesting 
an alternate model where they do not survive until today, but get converted into 
ellipticals by coUisions. 

Subject headings: galaxies: formation — galaxies: evolution — galaxies: high- 
redshift — galaxies: irregular 

1. Introduction 

A high fraction of young galaxies show an irregular clumpy structure from blue emission 
that is concentrated into clusters much larger and more massive than any found locally. 
Chain galaxies, for example (Cowie & Hu 1995; van den Bergh et al. 1996; Moustakas et al. 
2004) are dominated by several giant blue clumps strung along a line; there are no exponential 
light profiles or central red bulges as in spiral galaxies today (Elmcgreen, Elmegreen, & Sheets 
2004, hereafter Paper 1). Dalcanton & Schcctman (1996) suggested that chains arc edge-on 
low surface brightness galaxies. O'Neil et al. (2000) considered high-redshift disk galaxies 
in general, many of which have bulges and exponential profiles, and concluded from the 
distribution of axial ratios that chains are edge-on disks. Reshetnikov, Dettmar, & Combes 
(2003) showed that the axial ratios of 34 edge-on galaxies in the Hubble Deep Field are larger 
than for local galaxies, implying thicker disks, and noted that half of their sample, including 
some face-on galaxies, have non-exponential profiles like the chains. Edge-on local galaxies 
resemble chain galaxies too when seen in projection at ultraviolet wavelengths, although the 
local star forming regions are much smaller than in chain galaxies (Smith et al. 2001). 

Elmegreen, Elmegreen & Hirst (2004a, hereafter Paper II) studied a large sample of chain 
and non-exponential clumpy galaxies in a deep Advanced Camera for Surveys image and 
confirmed that they are members of a new class of young galaxies with thick, non-exponential 
disks dominated by ~ 10^ Mq blue clusters. They have no regular bulge structures, like other 
disk galaxies at the same redshifts (Paper I), but sometimes have irregular or offset bars with 
a bulge-hke size (Elmegreen, Elmegreen, & Hirst 2004b, hereafter Paper HI). An objective 
name for this type was suggested to be a "clump cluster," because without kinematic data 
it was not known whether each is a conglomeration of pre-existing small galaxies, like a 
tiny Stephan's Quintet, or a single coherent galaxy. A lack of clear rotation in some linear 
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galaxies at high redshift suggests they could be conglomerates (Bunker et al. 2000; Erb et 
al. 2004). Chain galaxies also lack rotation in the models by Taniguchi & Shioya (2001), 
where the objects are shocked filaments between expanding galactic-wind bubbles, and by 
lono, Yun, & Mihos (2004), where they are interaction- induced bars with gas inflow. If they 
are in fact coherent objects, then "clumpy galaxies" might be a better term. Here we refer 
to them as clump-clusters or clump-cluster galaxies, depending on the context. 

Conselice et al. (2004) noticed a similar dominance of clumps in many yoimg disk galax- 
ies, including galaxies with bulges and exponential-like proflles, and called them "luminous 
diffuse objects." Roche et al. (1997) found that faint and large blue galaxies with low 
surface brightness tend to have spiral or chain-like morphologies. Paper 1 showed that at I 
band, chains dominate all edge-on galaxy types by a factor of ~ 2 beyond 24th magnitude. 
Conselice et al. (2004) showed that luminous diffuse objects dominate star-forming galaxies 
at redshifts z > 1.5. Conselice, Blackburne, & Papovich (2004) found that peculiar galaxies 
dominate at z > 2. Consistent with this, there are very few nearby chain galaxies having 
I< 22 mag (Abraham et al. 1996). 

To study further the formation of stars and disks in young galaxies, we surveyed the 
Ultra Deep Field (UDF; Beckwith et al. 2004) for more examples of chains and clump- 
clusters. We identified 178 clump clusters with major axes exceeding 10 pixels, and measured 
their integrated colors, magnitudes, and surface brightnesses as well as the properties of the 
individual clumps and interclump regions. The lack of exponential disks was confirmed with 
contour plots and radial profiles. Visual inspection is preferred for clumpy galaxies because 
the SExtractor program sometimes mistakes them for multiple sources (Paper I). Here we 
discuss 10 examples of this class. A full catalog will be presented in Elmegreen et al. (2005). 
The UDF is processed for a resolution of 0.03 arcsec per pixel; at a redshift between 1 and 
3, 0.03 arcsec corresponds to a spatial scale between 220 pc and 230 pc in a ACDM model 
with Qa = 0.73, Qm = 0.27, and Hq = 71 km s"^ Mpc"^ (Spergel et al. 2003). 

Galaxy evolution models of the colors and magnitudes of distant galaxies have been 
used to estimate the redshifts, ages, extinctions, and stellar masses. Studies of galaxies in 
the Hubble Deep Fields North (Papovich, Dickinson Sz Ferguson 2001; Dickinson et al. 2003) 
and South (Rudnick et al. 2003), the Chandra Deep Field South (Idzi et al. 2004), and in 
fields observed from the ground (Shapley et al. 2001) show a general trend of decreasing star 
formation rates and increasing stellar mass over time to the present (sec review in Ferguson, 
Dickinson, & Williams 2000). We are interested in the ages and masses of the clump and 
interclump regions in our sample and use evolution models like these to compare with the 
observations. We would like to determine if the clumps formed in initially smooth disks 
as a result of spontaneous gaseous instabilities (Noguchi 1999; Immeh et al. 2004a,b), or 
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if they formed by some other mechanism, such as compression-triggered star formation in 
giant, hierarchically-merging clouds (Brook et al. 2004). The clumps could also be separate 
galaxies combined during minor mergers (Walker, Mihos, & Hernquist 1996; Abadi et al. 
2003). The instability model might apply if the gas accretion was relatively smooth, as in 
galaxy formation models by Murali et al. (2002), Westera et al. (2002), Sommer-Larsen, 
Gotz, & Portinari (2003), and Keres et al. (2004). Smooth gas accretion on the scale of the 
disk could occur in the hierarchical model if the gaseous parts of the accreting galaxies are 
stripped from the stars at large distances from the disk and then settle hydrodynamically. 

Consclice (2004) suggested that the early build-up phase of galaxy formation, corre- 
sponding to 2; > 3, is dominated by major mergers, the intermediate phase at 1 < ^ < 2 
is dominated by hydrodynamic gas accretion, and the late phase {z < 1) is dominated by 
minor mergers. He noted that Luminous Diffuse Objects, which are morphologically similar 
to clump clusters, have bright outer regions with low central light concentrations, and so 
might be examples of the intermediate phase. The signature for this may be present in the 
age and mass distributions of the clump and interclump regions, studied here. 

2. Data Analysis and Results 

The UDF images used for this study were observed and processed by Beckwith et al. 
(2004) and are available on the STScI archive. The UDF consists of images in 4 filters: 

F435W (hereafter B435 band; 134880 s exposure), F606W (Veoe band; 135320 s), F775W 
(1775 band; 347110 s), and F850LP (zgso band; 346620 s). The magnitude zero points are 
posted on the archival data sites: 25.673 (B), 26.486 (V), 25.654 (i), 24.862 (z). The images 
were mosaiced and drizzled by the HST team to produce final images of 10500 x 10500 pixels 
with a scale of 0.03 arcsec per px. Beckwith et al. used the program SExtractor to identify 
some 10,040 sources, so we use their nomenclature for the sources we found independently 
by visual examination. 

Figure 1 shows 10 of the most spectacular clump clusters seen at i775-band. They were 
chosen because of their large angular size, good resolution into clump and interclump regions, 
and wide range of feature types representing most galaxies in this class. Their positions and 
integrated properties are in Table 1. Some have multiple UDF numbers because SExtractor 
identified the individual large clumps as separate sources. Contour plots of intensity are 
drawn with the lowest contour la above sky noise, corresponding to a surface brightness of 
26.8 mag arcsec"^ at iyys-band. The lower left-hand bar in each panel represents 0.5 arcsec. 
The numbers in the lower right are the UDF catalog numbers. 
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Objects UDF 3483 and UDF 82+99+103 are ring-like galaxies, while UDF 7081+7136+7170, 

UDF 4245+4450+4466+4501+4595, UDF 3034+3129, UDF 3465+3597, and UDF 3483 have 
linear central sections resembling bars (Paper III). UDF 4245+ also has several small clumps 
distributed around the periphery, while galaxy CC-1 has many small clumps throughout the 
body. UDF 2012 apparently has a dust lane running through the middle. 

The integrated galaxy magnitude was measured in a rectangular box with edges defined 
by the 1775 isophotal contours 2a above the sky noise for each filter (corresponding to a 
surface brightness of 26.0 mag arcsec^^). We also did photometry on each prominent clump 
and several interclump regions, giving surface brightnesses and colors. Clump photometry 
was done using imstat in IRAF by measuring the total fiux inside a box that surrounds the 
clump and then subtracting the adjacent background in the same size box. The background 
was based on several fiux readings in the immediate neighborhood of the clump. 

Bright clumps have surface brightnesses that are ~ 1.5 mag arcsec^^ brighter than the 
interclump regions. On average, ~ 40% of the i-band luminosity is in the clumps, as indicated 
by the clump luminosity fractions, /lum, in Table 1. This luminosity fraction corresponds to 
a different rest-frame wavelength in each case because the galaxy redshifts differ, so it is only 
meant to be a indicator of the highly clumped morphology viewed at i-band. The two cases 
with small /lum, UDF 4245+ and CC-1, have many more clumps, and each is much smaller 
than in the other objects. Most clump colors within a given object arc similar to each other, 
although some apparently dusty interclump regions are much redder than others, such as 
the dust lane in galaxy UDF 2012. 

A color-magnitude diagram for 1775 versus 1775-2350 is shown in Figure 2 for all of the 
big clumps in Figure 1 (circles) and for the sum of these clumps (small squares) and the 
integrated galaxies (large squares). The 1775 magnitudes of the clumps range from 25.4 to 
30.5. The colors lie mostly in a narrow range for each galaxy but they vary from galaxy to 
galaxy with the redshift (see below). 

Color-color diagrams for each clump-cluster galaxy are in Figures 3 and 4. The circles 
are the clumps, the crosses are the interclump media at various positions, the big squares 
are the integrated galaxies, the small squares are average colors of the clumps (weighted by 
clump luminosity), and the small triangles are the average weighted colors of the measured 
interclump emissions. Superposed curves are best-fit models that will be described in Section 
3. The reddest B435-V6O6 colors are in UDF 3034+3129 and UDF 2012. In UDF 3034+3129, 
the central clump in the image has B435-V6O6 — 1-8 mag, in UDF 2012, most of the red 
interclumps at B-V> 1 mag are in the dust lane. Also very red (B435-V6O6 = 1-1 mag) is the 
large western clump in UDF 4245+, which has a similarly red interclump region. The two 
bluest clumps in B435-V6O6 are in UDF 7081+ (the southernmost round clump) and UDF 
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2012 (the central clump just above the dust lane). The bluest clump in ^775 — zgbo is in CC-1, 
and is a faint clump in the southeast. Figures 3 and 4 show that the interclump regions are 
systematically redder than the clumps. 

3. Comparison with models 

3.1. Model Components 

Photometric redshifts for the 10 clump-cluster galaxies and the ages and luminous 
masses of the clumps and interclump regions were determined by comparison with stellar 
population models. The luminous mass is defined to be the mass that fits the observed emis- 
sion in the bandpasses we have available. This mass does not include possible contributions 
from older populations of stars that do not emit significantly in the observed bandpasses and 
are not modelled explicitly in the star formation history. 

The models used the GALAXEV tables of spectra in Bruzual & Chariot (2003) for the 
Padoval994 evolutionary tracks with a Chabrier (2003) IMF (extending from O.IM0 to 100 
Mq with a flattening below ~ 0.5 Mq). These tables give integrated population spectra at 
1221 wavelengths from 91 A to 160 /im and at 221 times from 1.25 x 10^ years to 2 x 10^° 
years. We use the low-resolution version. The spectra were integrated over time with two 
basic time dependencies for star formation, a constant rate, and an exponentially decreasing 
rate with timescales of 10^ 3 x 10^ 10^ 3 x 10^ and 10^ years. A choice of six metallicities is 
available in GALAXEV. Shaplcy et al. (2004) found that the metallicities for 7 disk galaxies 
at the same redshifts as those found here (2; ~ 2 — 2.5) arc approximately solar {Z = 0.02). 
The metallicity of the Milky Way thick disk, wliich may have formed at about this epoch 
(Dalcanton & Bernstein 2002; Brook et al. 2004) peaks at [Fe/H] ~= -0.6 {Z = 0.005). We 
use Z — 0.008 in what follows; different metallicities did not change the results significantly. 

Absorption from cosmological hydrogen lines was included according to the prescription 
in Madau (1995), taking Lyman series lines up to order 20 and including continuum absorp- 
tion. Figure 5 shows the resulting transmission for a source at 2; = 3.5 for comparison with 
Figure 3 in Madau (1995). Dust absorption was included using the wavelength dependence 
in Calzetti et al. (2000) with the short-wavelength modification in Leitherer et al. (2002) 
and the redshift dependence for galaxy intrinsic Ay extinction in Rowan- Robinson (2003). 
The resulting rest-frame dust transmissions for galaxies at z = 1 and z = 3.5 are shown 
in Figure 6. Other estimates of extinction in 2; ~ 2 disk galaxies arc 2 or 3 times higher 
than these (e.g., Adclbcrgcr & Steidel 2000; Papovich, Dickinson, & Ferguson 2001), so two 
other sets of models with twice and four times the Rowan-Robinson extinction were fit to 
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our data for comparison, as discussed below. In general, higher extinction lowers the age 
and mass-to-luminosity ratio as the intrinsic stellar populations become bluer and brighter, 
but the masses may either increase or decrease depending on the relative size of the changes 
in M/L and 10*^-^"^ for extinction A; the redshift solutions do not change noticeably with 
extinction over this range. 



3.2. Model color evolution 

The time-integrated spectra, corrected for cosmological hydrogen absorption and intrin- 
sic galactic dust, were convolved with the four throughput functions for B435, Vgoe, 1775, and 
Zgso filters, as determined from the IRAF task SYNPHOT (kindly run for us by Jennifer 
Mack at STScI). Equation 8 in the postscript file supplement of the GALAXEV download 
was used in the AB magnitudes convention. Because this AB version was not written ex- 
plicitly for intensity as function of wavelength, we give it here: 



m — —2.5 log 



48.60. (1) 



To get the constants of proportionahty, the luminosity L in GALAXEV, which is tabulated 
in units of Lq for each solar mass of stars, was converted to erg s~^ using Lq = 
3.826 X 10^^ ergs s~^; c is the speed of light in cm s~^, di is the luminosity distance in cm. 
The filter function, -F(A), is dimensionless, as is d\/\ in the denominator. The combination 
d\L\ has dimensions of ergs when A is in A, but the other A in the numerator has to be 
converted to cm to give A/c = 1/ (3 x 10^^) s^; A is the observed wavelength. 

Models of color versus time are in Figure 7 for three values of redshift and for the 
six star formation histories. The lower curve in each panel (bluest color) is for continuous 
star formation and the upper curve (reddest color) is for the most rapidly decaying star 
formation, with an exponential decay time of 10 My. Intermediate curves span the range 
of exponential decay times, r, indicated in the figures. The colors rise with time rapidly, 
especially at high redshift, because of the depletion of uv hght in the absence of current star 
formation. The Lyman absorption jump shown in Figure 5 reddens the B435-V6O6 values for 
all of the curves when z > 3.5 (the drop-out effect - Steidel et al. 1996). The redshifts of 
star formation corresponding to the times on the abscissae are indicated on the top axes. 
Analogous models without intrinsic dust extinction differed only a little in color, at most 
~ 0.1 mag in B435-V606- 



Model color-color plots are shown in Figures 8 and 9. Stellar age in the rest frame of 
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the galaxy increases along the curves from at the bottom left (bluest colors) to whatever 
age corresponds to z ~ 8 at the observed redshift (5.3 Gy, 2.7 Gy, and 1.5 Gy at ^ = 1, 
2, and 3, respectively). The different star formation histories have separate curves. Solid 
curves have dust extinction from Rowan-Robinson (2003) and dashed curves have 4 times 
this extinction. Generally, the more continuous the star formation, the closer the curve stays 
to the lower left end, even at large times. The most extreme red values occur for the most 
rapid decay, where the exponential time is r = 10 My. 



3.3. Photometric redshifts 

The observed colors of the clump and interclump emissions correspond to the lower 
left parts of the theoretical color-color curves. Figures 3 and 4 show the data and models 
together, using the redshift models that best match the colors. To find these redshifts, 
we made models with the six star formation histories mentioned above and with redshifts 
ranging from z — to 6.75 in intervals oi z — 0.25. Model colors were tabulated for a wide 
range of times within each history type. The rms deviations between the average clump 
positions and the models on both color-color plots were then determined for each clump- 
cluster galaxy, and the model with the minimum rms was chosen. The highly deviant colors 
shown by arrows in the figures were not used for this. Thus the photometric redshifts are 
based only on the clump colors, not the average galaxy colors or the interclump emission. We 
do this because clumps are more homogeneous than whole galaxies. The fit to the redshift 
also gives the preferred star formation decay time, r, and the average age of the clumps in 
each object. 

We solved for the three main variables (redshift, decay time, age) in a certain way 
to minimize possible variations from intrinsic scatter in the data. We first solved for the 
minimum value oi 5 — {bv — hvmY + {vi — vimY + iiz — iZmY along each curve that traces out 
age for a given redshift and star formation history. This minimum was found by a quadratic 
fit to the lowest three 5 values along that curve. Here, hv, vi, and iz represent the observed 
colors and hvmi vim and izm represent the model colors. Next, the minimum values of 5 were 
arranged in order of increasing redshift and, for each redshift, increasing log(r) in the star 
formation history. Along this list there were usually one or two broad minima in the 5 values, 
corresponding to close fits between the observed and model colors. Typically one close fit 
was at low z and another was at intcrmediate-to-high z. This duplicity corresponds to the 
looping structure of the color-color curves in Figure 8, i.e., the template spectra have the 
same slope at two distinct wavelengths (typically at the rest optical and rest uv) for some 
star formation histories and extinctions. We could rule out the small-2; solutions because 
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these gave effective rotation curve speeds (see below) tfiat are unreasonably low for disk-like 
galaxies (i.e., 20 km s~^ or less). The broad minima along the ordered models were then 
isolated from each other by considering only S values less than certain limits, 0.1, 0.05, and 
0.02. For each limit, the mean and the standard deviation of the z value was determined, 
weighted by the inverse of 6. When this mean z was about the same for each limit, and 
the standard deviation in z became small, ±0.2 or less, for the smallest limit, then we felt 
confident in our z solutions. 

The redshifts based on the whole galaxy colors were typically several tenths smaller than 
the redshifts based on the average clump colors. This difference is larger than the standard 
deviation for each z and also larger than the differences in fitted z's for the three limiting 
5 values. Most likely, the interclump emission, which contains about half the total light, is 
contributing to the composite galaxy spectrum in a way that makes the composite colors 
inconsistent with a single exponential-decay model for star formation. For example, whole 
galaxies are redder than clumps in irrs-zssO) and the models in Figure 8 shift 1775-2350 toward 
the red for lower z. We believe that the fits to the clump colors give more accurate redshifts 
as long as we use simple star formation histories. 

For each fitted z, we also got the inverse-(5 weighted mean of log(T) in the star formation 
history and the inverse-5 weighted mean of the age of the star formation region, t. 

To check our method, we compared the photometric redshifts derived from our models 
with the 19 pubhshed redshifts available for the GOODS field (Vanzella et al. 2004) using 
the photometry in that study and the same filters as the UDF field. We also compared 
our redshifts to those derived by Benitez et al. (2004) in the Tadpole Galaxy field, using 
the different filters appropriate to the Tadpole Galaxy survey. Figure 10 shows a plot of 
these redshift distributions. Our values matched the GOODS redshifts, which range from 
0.21 to 1.4, with an rms deviation of 0.2, not counting an object at z = 5.83, for which we 
derived z = 2.5. This single large discrepancy could result from strong emission lines in the 
observation that are not present in the models, or from inappropriate models at high z. Our 
values matched the 107 high-quality photometric redshifts in the Tadpole field (i.e., those 
for which ODDS> 0.99) to within an rms of 0.3, not counting 29 more for which Benitez et 
al. derived a redshift of around 1.85. These excluded galaxies (not plotted) form a cluster in 
the redshift distribution with a value of z from Benitez et al. ~ 1.872 ± 0.015; our redshifts 
for them range from 0.3 to 2.1, with an average of 1.28 ±0.43. Aside from these unexplained 
discrepancies, our photometric redshifts are comparable to the redshifts derived by others 
to within several tenths. We also determined redshifts for 30 B-band dropout galaxies in 
the UDF field (Elmegreen et al. 2005) and obtained an average of 3.49 ± 0.37, which is 
reasonable. 
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The interclump age could not be derived reliably from our models because of the am- 
biguity between age and decay time r. The interclump emission is probably a composite 
of many ages, however, going back to the beginning of Pop II star formation, regardless of 
the nature of the galaxies in which these stars formed. To find the best interclump star 
formation histories, we used the galaxy redshifts found from the clump colors and assumed 
an interclump age corresponding to a beginning of interclump star formation ai z = 6. We 
then determined the value of r for each galaxy that best fits the average three interclump 
colors. 



3.4. Redshift and age solutions 

Table 2 shows the main results of the fitting procedure with the Rowan-Robinson ex- 
tinctions. The average clump age is derived to be 0.34 Gy from the clump colors and models. 
The average interclump age is assumed to be the time between the Universe age at the fitted 
clump redshift and the Universe age at 2; = 6, when the interclump star formation presum- 
ably began. The average value in Table 2 is 2.14 Gy. This larger age for the interclump 
region is not unreasonable considering the redder colors there. The clumps are characterized 
by a decaying star formation rate with an e- folding time that has an average r ~ 0.1 Gy. 
Thus the clumps are three times older than their e-folding times and should not be con- 
sidered starbursts. The fitted e-folding times for star formation in the interclump regions 
average ~ 1 Gy, which is only slightly less than the average interclump age and consistent 
with a continued build-up of these regions. 

Extinctions E(B-V) larger than the Rowan-Robinson (2003) values by factors of 2 and 4 
for both clump and interclump emission give about the same redshifts for the galaxies. The 
average values for various quantities in these two cases are given in the bottom rows of the 
tables. The average redshift decreases from 2; = 2.3 with our fiducial extinction to 2.2 for 
both 2x and 4x extinctions. The higher extinctions decrease the average clump ages from 
0.34 Gy to 0.23 Gy and 0.09 Gy, respectively, with < r >= 0.17 and 0.16 for the two cases, 
because the source has to be intrinsically bluer when there is additional reddening. They 
also make the interclump regions intrinsically bluer by increasing the star formation decay 
times, T, from an average of 1.01 Gy to averages of 2.09 Gy and 3.30 Gy, respectively. 

The redshifts allow the determination of the luminosity distance dL in equation 1 (Car- 
roll, Press, & Turner 1992). The resulting model apparent magnitudes scaled to the observed 
apparent magnitudes give the clump masses. This derivation of mass uses the evolution 
models in GALAXEV again, as they give the residual stellar mass per initial solar mass 
of a population of stars as a function of time. To get the clump mass, we integrated this 
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residual stellar mass over time using the same time-dependent star formation rate as for the 
integrated spectrum, and then multiplied the resultant model mass by the factor used to 
scale the model luminosity to the observed clump luminosity 

3.5. Interclump mass column densities 

The interclump mass column densities were obtained from the observed interclump 
surface brightnesses. To be careful about surface brightness dimming, we first wrote the 
absolute magnitude of a solar mass of stars from equation 1 using di = 10 pc. Wc then 
used the cosmological distance modulus (Carroll, Press, & Turner 1992; Paper 111) to give 
the apparent magnitude for this solar mass and scaled this result to the observed apparent 
magnitude in a square arcscc. This gave the mass per square arcscc. The square arcsec 
was then converted into a square kiloparsec using the cosmological angle formula (Carroll, 
Press, & Turner 1992) to get the rest-frame mass column density. This process corrects 
for the z) surface brightness dimming. We got the same result using the z-dependent 
luminosity distance dz, in equation 1. The total interclump mass was determined by averaging 
the interclump masses per square arcsec for each galaxy and multiplying this average by 
the galaxy area in square arcsecs. This assumes the interclump medium is uniform. The 
interclump mass column density is typically larger than the clump mass column density. 
Evidently, the clump light is dominated by the brightest blue stars, while the total mass in 
the region is dominated by the underlying old stars. The sum of the clump mass and the 
interclump mass is taken to be the luminous galaxy mass. 

3.6. Results 

Table 3 gives the average clump mass, < Mdump >) diameter < -Dciump >■, and clump 
mass fraction, /mass, in each clump-cluster galaxy, and it gives the luminous galaxy mass, 
Mgai, the galaxy diameter, D^si (equal to the major axis length out to the 2-a contour), and 
the effective rotation speed, v^ot = [GMgai/ (5L'gai/2)]^'^^ (the factor 5 in the denominator 
assumes a uniform spherical mass distribution). For the physical sizes, D, the angular sizes 
were used along with the cosmological angle-size relation for the galaxy redshift. 

The results in Table 3 may be summarized as follows. The objects we have selected are 
large composite systems with an average luminous mass of 6.5 x 10^° Mq and an average 
diameter of 19.1 kpc out to the 2a contour. The size and mass are not unprecedented for 
disk galaxies at redshifts of 2 to 3 (Labbe et al. 2003; Shapley et al 2004). The average 
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rotation speed would be 150 km s"^ if there were regular circular motion, but there are no 
spectroscopic observations yet. There is a range of a factor of 50 in mass and a factor of 8 
in rotation speed, so the peculiar morphology of this galaxy type does not correspond to a 
particular mass. The clumps represent an average total mass fraction of 19%, again with a 
factor of 10 variation. This is smaller than their luminosity fraction mentioned above (40%) 
because the older interclump medium has a higher mass-to-light ratio than the clumps. 

The average interclump mass column density is ~ 4 times larger than the mass column 
density in solar neighborhood. The primary reason for the high value is the red color of the 
interclump stars. Colors this red are matched to stellar populations with high mass-to-light 
ratios. Increasing the assumed dust abundance by factors of 2 and 4 for both the clumps and 
the interclump regions hardly changed the redshifts and also changed the clump masses only 
a little, from an average of 0.64 x lO'' Mq to averages of 0.97 x lO'' Mq and 0.53 x 10*^ Mq, 
respectively. These small changes reflect the close balance between decreasing M/L ratios 
for intrinsically bluer stars and increasing exp (raust) corrections for more reddened stars. 
The interclump masses and total galaxy masses changed with extinction only a little too: 
the average total galaxy masses changed from 6.5 x 10^° Mq to 6.4 x 10^° Mq 11 x 10^° M0, 
respectively. Lower metallicity would also make the intrinsic colors bluer, and this would 
make the interclump ages older to match the observed colors, producing an even higher 
M/ L. Changes in the stellar initial mass function could lower the inferred interclump mass, 
but only if the relative abundance of stars younger than ~ 1 Gy, corresponding to A-type 
stars or earlier, is higher in proportion to low mass stars than found locally; i.e., the IMF 
has to be relatively flat at intermediate to high mass. 

In another series of models, we assumed the Rowan-Robinson (2003) extinction for the 
models of the clump emission, getting the same clump redshifts and ages as before, but 
then we assumed higher extinctions in the interclump regions by factors of 2 and 4 before 
deriving their star formation histories and masses. These higher interclump extinctions gave 
intrinsically bluer colors and a greater extinction correction to the magnitudes. The resultant 
interclump masses averaged for all 10 galaxies increased from 5.9 x 10^° Mq to 6.5 x 10^° Mq 
12 X lO^*' Mq respectively, as the corresponding lower mass-to-luminosity ratios were more 
than offset by a higher mass needed to overcome the greater extinction. 

Other clump properties are in Table 4. The clump equivalent atomic density, ndump, is 
taken to be the ratio of the clump stellar mass to the volume given by the clump diameter, 
-Dciump; assuming a mean molecular weight 2.2 x lO^^^g. This equivalent density is useful for 
comparisons with molecular cloud and other gas densities (see below). The dynamical time 
is tdyn = {Gp) for mass density p. The average and current star formation rates (SFR) 
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are given by 

SFR(ave) = M^.^^Jt ; SFR = (2) 

for clump age t and decay time r in Table 2. The internal clump velocity dispersion assumes 
virial equilibrium with a uniform density: Afdump = [G^^ciump/ (2.5Dciump)]^''^- The last 
column is the ratio of the clump internal dynamical time to the galaxy dynamical time. All 
of these properties are discussed below. 



4. Discussion 

4.1. General considerations 

Clump-cluster galaxies are a new type of object seen primarily in deep fields such as 
those surveyed by HST. They are characterized by: (1) a disk-hke shape (Paper II); (2) a 
significant contribution to the total luminosity (~ 40% - Table 1) from a small number of 
relatively blue clumps; (3) a fiat hght profile along the major axis (Papers I, II); (4) no central 
red clump or bulge (Papers I, II); and (5) a faint, relatively red interclump emission (this 
paper). Some clump-clusters also have bar-like features (Paper III). Edge-on clump-clusters 
may be chain galaxies (Paper II). The lack of clear rotation in chain galaxies (Bunker et al. 
2000; Erb et al. 2004) is a mystery if these are normal galaxy disks. 

In the few cases studied here, the clumps are typically ~ 0.3 Gy old in their rest frame 
and not currently forming stars as fast as they once did. Emission from the interclump 
medium is dominated by stars that are much older, perhaps up to 2 Gy older in some cases. 
Older and fainter populations inside the clumps could be present but undetected. 

With ~ 19% of the himinous mass in a dozen or so giant clumps, these galaxies are unlike 
modern galaxies where the largest star forming regions are typically < 10^^ of a galaxy's 
mass. This difference implies several things. First, the extreme dumpiness suggests that the 
rotation curve will not be simple (Immeli et al. 2004a,b). The clumps should attract each 
other and scatter, with their coldest and most dissipative parts migrating to the center to 
make a bulge and their hottest, least dissipative parts filling the interclump region or thick 
disk (Noguchi 1999; Somerville, Primack & Faber 2001; Immeh et al. 2004a,b; Hartwick 
2004; Brook et al. 2004). The lack of clear rotation in the few chain galaxies that have been 
observed may be the result of this severe clump interaction. 

Second, the clump ages are longer than their dynamical times, so the clumps are self- 
bound like giant clusters. The clump density is comparable to the mean galaxy density, 
however, so tidal forces on the clumps are significant (this follows from the comparable clump 
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and galaxy dynamical times in Table 4). This means that galactic shear and tidal stretching 
are moderately important for the lowest density clumps. For example, UDF 6462+6886 and 
UDF 82+99+103 have clumps with tails and arcs suggesting the influence of tidal forces. 

The fact that the equivalent clump atomic densities are much less than typical molecular 
cloud densities implies that the clump stars probably formed in numerous dense cores inside 
kpc-size, gravitationally-bound, low-density cloud complexes. This hierarchical morphology 
for star formation is analogous to the case of star complexes found in local spiral galaxies 
(Efremov 1995), except that for the clump clusters, the clumps are more strongly bound; 
local star complexes readily disperse or shear into flocculent spiral arms. 

The internal velocity dispersions in the clumps are inferred to be moderately large, 24 
km s~^ on average, based on the virial theorem. This velocity is ~ 16% of the average 
galactic orbital speed - a high fraction for gas in local spirals but comparable to the fraction 
in dwarf galaxies. If the clouds that made the young stellar clumps formed by gravitational 
instabilities (Noguchi 1999; Immeh et al. 2004a,b) and if they still have the ambient Jeans 
mass, then the velocity dispersion in the pre-cloud gas had to be comparable to the virial 
velocity in the clumps (this is based on the usual isothermal assumption). Such a high 
dispersion implies that clump-cluster galaxies are far more turbulent than local galaxies. 

One source of turbulence could be major encounters with other galaxies. Tidal forces 
increase the velocity dispersion considerably in local galaxies undergoing such encounters 
(Kaufman et al. 1997). The galaxies in our sample have no major companions now, but 
the giant clumps in them could have been triggered ~ 0.3 Gy ago at previous encounters, 
considering the average clump age. Another possible source of turbulence is accretion from 
the intergalactic medium. The disk impact speed can be ~ 0.4 times the circular speed if 
the gas falls in along cosmological filaments (Keres et al. 2004). This would be ~ 60 km s^^ 
for our galaxies. Disk turbulence should be slower than this because of energy dissipation, 
but still enough to explain the high Avdump- Supernovae could also be an important source 
of turbulence as they are in local galaxies, but in clump-cluster galaxies they are likely to 
be highly concentrated in the clumps, and then they may just cause blow-out rather than 
widespread turbulence. 

In the accretion model, clump-cluster galaxies would have added intergalactic gas at 
a rate comparable to the average star formation rate, ~ 20 yr~^. The disks would 
presumably have become unstable when the mass column density exceeded a critical value. 
Considering the clump ages, the resulting star formation has a timescale of ~ 0.3 Gy, which 
is several galactic orbit times. 
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4.2. Accretion models 

A different possibility is that the clumps were formerly separate galaxies that accreted 
during hierarchical buildup (e.g., Abadi et al. 2003; Brook et al. 2004). Some of our clump 
clusters are so irregular in morphology (e.g., UDF 6462+6886) that they appear more like 
random conglomerations than circular disks. Also, the lack of clear rotation in the few 
cases observed suggest highly disturbed dynamics. This model raises three questions: (1) 
Can small galaxies like our clumps survive their capture by large galaxies, considering the 
interclump population that was present before? (2) Are there examples of bare chimps in 
the UDF that resemble the clumps in our sample? (3) Would the distribution of axial ratios 
for coalesced systems be the same as that observed for clump clusters and chain galaxies? 
We attempt to answer these questions here. 

4-2.1. Clump Survival 

The short-term survival of a small galaxy captured by a large galaxy depends mostly on 
the ratio of densities. This ratio is important because it is also the ratio of the self-binding 
force per unit length inside the small galaxy to the tidal disruptive force per unit length 
exerted by the large galaxy. If the small galaxy is much denser than the large galaxy, then 
the small galaxy will survive until it migrates, perhaps because of dynamical friction, to a 
sufficiently dense region of the large galaxy, at which time it disperses. The migration time 
depends on the orbit, and the outer parts of the clump can be shredded continuously as it 
migrates. In the clump-cluster galaxies studied here, the typical column density of a clump is 
comparable to the average column density of the interclump population, both being several 
hundred Mq pc~^. The clump density should be considerably higher than the interclump 
density, however, because the interclump medium is likely to be much thicker than any single 
clump as a result of stirring from clump motions, dissolution of previous clumps, and stellar 
migration. The ratios of axes for these galaxies also suggest fairly thick disks (Paper II, and 
see below). In addition, the clumps probably have dense cores like other bound clusters. 
The clumps are therefore likely to survive their first few impacts with the disk, but probably 
not much longer. The observation that the clumps are much older than both their internal 
dynamical times and the galaxy orbit times (by a factor of 10) implies that most of them live 
outside the densest regions of the disk. The exceptions are the clumps with morphologies 
suggesting tidal disruption, such as the head-tail shapes of the clumps in UDF 82-1-99-1-103, 
UDF 6462+6886, and UDF 3034+3129. 
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4-2.2. Bare clumps in the UDF 

What is the evidence for clumps hke these outside of clump-cluster galaxies, where 
they presumably existed before the accretion? The UDF contains many isolated objects 
that resemble our clumps in both luminosity and color. The right-hand side of Figure 11 
shows color-magnitude and color-color plots of all the UDF objects smaller than 6 pixels 
in FWHM, as given in the tabulation on the UDF web site^. This size limit was chosen 
because it roughly corresponds to the size of a clump in a clump-cluster galaxy. The left side 
of the figure shows the magnitudes and colors for the measured clumps in the clump-cluster 
galaxies of this paper. The distributions are essentially the same in the regions of overlap. 
The UDF clumps can be much bluer than our clumps, suggesting either that star formation 
begins to slow down once a clump is ingested into a larger galaxy, or the clump-cluster 
sample in our survey has too few clumps to include the rare active ones seen in the general 
UDF field. There are also much fainter clumps in the UDF field than those measured in 
our survey, but the clump-cluster galaxies have much fainter clumps too which we did not 
study. We note that these small objects in the UDF field are not the well-studied Lyman 
Break galaxies, which tend to be more massive, blue, and luminous than our clumps even 
at the same z ~ 2.5 — 3 (e.g., Papovich, Dickinson, & Ferguson 2001). The small field 
objects are more the size of the 2; ~ 6 galaxies studied by Bunker et al. (2004). Most likely, 
most are low-mass and low-luminosity galaxies in about the same range of redshifts as the 
clump-cluster galaxies studied here. 

4.2.3. Distribution of axial ratios 

The distribution of axial ratios for a random assembly of field clumps might be expected 
to be significantly different than the distribution for a thick disk. This distribution was 
determined for a combined sample of chain and clump-cluster galaxies in Paper II and 
shown to be similar to the distribution for nearby spiral galaxies (see also O'Neil et al. 
2000; Reshetnikov, Dettmar, & Combes 2003). In general, a thin disk should have a fiat 
distribution for the ratio of minor to major axes, W/D, with a slight rise and then a drop at 
low ratios where the projected width W becomes comparable to the intrinsic disk thickness. 
A galaxy whose appearance is dominated by a few clumps will not be very different, in fact. 
The first three clumps make a plane regardless of their positions, so the question is how 
much do the few additional clumps in a clump-cluster or chain galaxy distort this projected 
plane? 



^http: / /archive. stsci.edu/pub/hlsp/udf/acs-wfc/h_udf_wfc_ VI_I_cat.txt 
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Figure 12 shows distribution functions of axial ratios for model galaxies with various 
numbers of clumps randomly positioned in 3 dimensions. Each line is a histogram from 10^ 
models. The first three clumps in each model are taken to define a plane and the others are 
given random positions around the plane in polar coordinates. The radial distance of each 
clump is distributed as a Gaussian. To simulate variations from different degrees of flattening 
in a disk, the distance perpendicular to the plane is multiplied by factors of Z = 0.1, 0.3, 
0.5, 0.7, and 1 for the five different lines. Once a clump distribution is chosen, the polar axis 
of the model is tilted by a random angle whose cosine has a uniform distribution function, 
simulating random inclinations of the model in the sky. The axial ratio used for the histogram 
is then taken to be the minimum ratio of projected width to length, considering all angles 
in the sky plane. The resulting histograms peak for axial ratios near unity when Z = 1 and 
there are a lot of clumps, because then the model galaxy is a uniformly filled sphere. When 
the number of clumps is small, the distribution decreases at high axial ratios because the 
clumps arc irregularly positioned even in the spherical case {Z = 1) and some directions 
happen to have have short sizes. The distributions decrease at small axial ratio because 
there are usually some clumps at non-zero perpendicular distance, making W/D non-zero. 
The axial ratio where this decrease begins is smaller for more fiattened systems because the 
minimum value of W/ D is ~ Z. 

The observed distribution of axial ratios for all of the chains and clump-cluster galaxies 
in our survey of the UDF field (Elmegreen et al. 2005) is also shown in Figure 12. The solid 
fine is for chains and the dashed fine is for clump-clusters, added on top of the solid fine. The 
observed distribution resembles the models only when there is some intrinsic fiattening of the 
galaxies. However, the degree of flattening does not have to be large. The best models are 
for a number of clumps appropriate to our galaxies, ~ 5 — 10, and for a flattening Z ~ 0.3, 
which is the second curve from the left (red in the electronic version). The corresponding 
intrinsic axial ratio of the model galaxies is the same, ~ 0.3, which implies the disks are very 
thick. The same conclusion was reached in Paper 11. 



4.3. Assessment of accretion model 

These considerations suggest that the clumps in some of our clump-cluster galaxies could 
have been accreted as whole objects from the field. They would probably survive for the 
required time, there arc isolated field galaxies like them, and the clump-clusters arc relatively 
thick disks. Some energy dissipation perpendicular to the disk had to occur to fiatten the 
clump distribution to the observed axial ratio of 0.1 to 0.3, but dynamical friction from the 
observed interclump disk could have done this during the capture process. 
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This capture scenario is consistent with the model by Brook et al. (2004), who suggest 
that thick disks in spiral galaxies form quickly, in less than 1 Gy, as a result of mergers with 
a small number of ~ 2 x 10^ Mq clumps. Hartwick (2004) also notes the quick formation 
time of thick disks, and suggests that the pre-galactic clumps which make them could remain 
today, in part, as disk globular clusters. 

A peculiarity of our sample is the high interclump surface density. If this result is 
correct, then it implies that the galaxies cannot simply age to form present-day spiral disks. 
Today's disks either have high mass column densities along with massive bulges, or low 
mass column densities without prominent bulges. They do not have high mass column 
densities and no bulges, like our clump-clusters. One possibility is that a large fraction of 
the mass in clump-cluster disks will migrate inward to form a bulge. This migration would be 
accompanied by outward angular momentum transfer, and might result from spiral waves or 
repeated interactions, with bar formation a possible intermediate step. Another possibility 
is that clump-cluster galaxies merge to form elliptical galaxies. The likelihood of this is high 
considering their relatively large sizes and high redshifts where interactions are still frequent. 

B.G.E. is supported by the National Science Foundation through grant AST-0205097. 
Support for D.M.E. was provided by the Salmon Fund of Vassar College. 
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Table 1. Properties of Galaxies 



UDF ID Numbers 


Position'' 


1775 


B436-V6O6 


V606-i775 


1775-2850 


A'775 


No. 


/lum 




3^ 32 


m 


-27° 


mag 








mag arcsec"^ 






82+99+103 


40.7472=, 


49' 


26.047" 


23.6 


0.43 


0.25 


0.30 


25.1 


8 


0.48 


3034+3129 


37.8404=, 


47' 


56.390" 


24.0 


0.82 


0.18 


0.08 


24.6 


5 


0.44 


2012 


35.9178=, 


48' 


17.216" 


24.3 


0.80 


0.11 


0.12 


25.0 


10 


0.38 


3483 


36.3842=, 


47' 


47.125" 


23.4 


0.20 


0.14 


0.13 


24.5 


12 


0.30 


3465+3597 


35.6573=, 


47' 


48.567" 


23.5 


0.11 


0.04 


0.02 


24.3 


11 


0.34 


CC-l'' 


49.3154=, 


47' 


26.377" 


22.0 


0.44 


0.16 


-0.07 


23.9 


11 


0.13 


4245+4450+ 


28.9470=, 


47' 


37.730" 


22.6 


0.08 


2.75 


-2.18 


23.7 


14 


0.06 


4466+4501+4595 






















1666 


41.3128=, 


48' 


21.154" 


23.9 


0.09 


0.10 


0.45 


25.2 


13 


0.39 


7081+7136+7170 


38.9909=, 


46' 


56.265" 


24.5 


0.06 


0.10 


0.34 


25.9 


8 


0.38 


6462+6886 


43.5596=, 


46' 


59.376" 


23.5 


0.18 


0.32 


0.27 


25.5 


8 


0.36 



^only the seconds of Right Ascension and the arc minutes and seconds of Declination are tabulated. 

''no UDF number given for this because it is near the edge of the image, so we use CC for "Clump Cluster". 



Table 2. Redshift and Star Formation Properties 



Galaxy 


redshift 


Clump 


Clump 


Intcrclump 


Interclump 


Number 


z» 


age, t (Gy) 


decay t (Gy)'' 


age (Gy)'= 


decay t (Gy) 


82+ 


2.5 


0.36 


0.11 


1.71 


0.33 


3034+ 


3.0 


0.21 


0.10 


1.24 


0.58 


2012 


3.0 


0.12 


0.03 


1.26 


0.65 


3483 


2.2 


0.26 


0.05 


2.02 


1.63 


3465+ 


2.4 


0.22 


0.06 


1.80 


1.06 


CC-1 


2.8 


0.35 


0.14 


1.40 


1.38 


4245+ 


1.7 


0.32 


0.05 


2.92 


1.16 


1666 


1.6 


0.81 


0.28 


3.15 


0.91 


7081+ 


1.6 


0.43 


0.10 


3.13 


0.92 


6462+ 


1.8 


0.31 


0.06 


2.82 


1.49 


averages = 


2.3 


0.34 


0.10 


2.14 


1.01 


ave, Av x 2 = 


2.2 


0.23 


0.17 


2.40 


2.09 


ave, Av x 4 = 


2.2 


0.09 


0.16 


2.32 


3.30 



''Photometric redshifts for the average clump colors 

''Average star formation history exponential decay rate for clumps 

■^Assumes interclump star formation began at z = 6. 
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Table 3. Derived Properties 



Galaxy 


< A^clump > 


^ ^clump ^ 


/mass 


A^gal 




frot 


Number 


109 Mq 


kpc 




lOlO Mq 


kpc 


km 


82+ 


1.36 


2.5 


0.04 


25 


18 


350 


3034+ 


0.70 


2.1 


0.16 


2.2 


16 


110 


2012 


0.17 


1.5 


0.07 


2.5 


13 


130 


3483 


0.75 


1.4 


0.32 


2.8 


14 


130 


3465+ 


0.50 


1.5 


0.28 


2.0 


17 


100 


CC-1 


0.96 


1.3 


0.07 


15 


25 


230 


4245+ 


0.21 


2.1 


0.59 


0.5 


22 


44 


1666 


0.52 


1.6 
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Table 4. Other Clump Properties 
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Fig. 1. — Grayscale i-band images and contour plots for 10 clump-cluster galaxies. The bar 
indicates 0.5"; galaxy UDF names are noted. The contours are linear in steps of 1 a sky 
noise, with the outermost contour at 26.8 mag arcsec"^. 
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Fig. 2. — Color- magnitude distribution for clumps (circles), for the sum of all clumps (small 
squares) and for the whole galaxy (big squares). Abbreviated galaxy names are indicated. 
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Fig. 3. — Color-color distribution of clumps (circles), interclump regions (crosses), and whole 
galaxies (squares). The average colors for the clumps (small squares) and interclump regions 
(triangles) are also shown. Curves are theoretical models that best match the fitted redshifts 
of the galaxies to the nearest 0.25. In general, the lower curves (bluer) correspond to the 
largest star formation decay times, which are assumed to be infinite (green curve for the 
electronic version of this publication), 10^ yr, 3 x 10^ yr, 10* yr, 3 x 10* yr, and 10^ yr 
(red curve). The best-fit models are determined simultaneously from both this color-color 
plot and the one in the next figure. Observed points that are far off the range of scales are 
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Fig. 4. — Color-color distribution of clumps (circles), interclump regions (crosses), whole 
galaxies (squares) and averages (small squares, triangles) are shown as in Figure 3, along 
with theoretical tracks that are closest in redshift to the fitted values. 
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Fig. 5. — Transmission from foreground cosmological hydrogen clouds for a galaxy at redshift 
z — 3.5 calculated using the method of Madau (1995). 




Fig. 6. — Transmission from dust inside the model galaxies using the Calzetti et al. (2000) 
and Leitherer et al. (2002) extinction curves and the intrinsic Ay extinction in Rowan- 
Robinson (2003). Dust produces less spectral structure and color variation than cosmological 
hydrogen absorption, so the redshifts derived here are not sensitive to internal extinction. 
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Fig. 7. — Model colors versus population age for three galaxy redshifts, z = 1, 2, and 3, 
and six models of the star formation history, as in Figure 3. Redshifts corresponding to 
the population ages for given galaxy redshifts are on the top axis. The ticmarks for these 
redshifts are not the same interval throughout, but they should still be self-evident; for 
example, on the top left, the ticmarks are at 1.002, 1.003, 1.004, etc., to 1.01, and then 1.01, 
1.02, 1.03, etc. to 1.1, and then 1.1, 1.2, 1.3, etc. to 2. 
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Fig. 8. — Population models on a color-color plot where time increases along each curve from 
the lower left to the upper right. The same six star formation histories as in the previous 
figures are used, with the same color coding. The solid lines are for internal extinction 
following the Rowan-Robinson (2003) result, and the dashed lines are for internal extinctions 
that are 4 times higher. 
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Fig. 9. — Model color-color plots as in Figure 8. 
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Fig. 10. — A comparison of rcdshifts calculated by different methods. On the bottom 
are photometric redshifts for the Tadpole Field galaxies calculated here (abscissa) versus 
the Benitez et al. (2004) photometric redshifts (ordinate). The filters appropriate for the 
Tadpole survey were used; 29 galaxies with Benitez et al. redshifts in a narrow range around 
1.8 were omitted as they appear to be anomalous. On the top are photometric redshifts 
for galaxies in the GOODS field (abscissa) versus the Vanzella et al. (2004) spectroscopic 
redshifts (ordinate). The filters appropriate to the GOODS survey were used in this case; 
they are the same as the filters for the UDF survey. 
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Fig. 11. — (left) Color-color and color-magnitude distributions for all the measured clumps 
in our 10 clump-cluster galaxies, and (right) analogous distributions for all the small clumps 
in the UDF field, with diameters less than 6 pixels. The similar distribution centroids suggest 
that the clumps in our clump-cluster galaxies could have been captured from the field. The 
faintest objects in the field also have faint counterparts in clump-cluster galaxies, but they 
were not included in our clump measurements. 



-34- 




0.2 0.4 0.6 0.8 

Axial Ratio, W/D 



Fig. 12. — (top) The distribution of the ratio of minor to major axes for all the chain galaxies 
(solid line) and clump-cluster galaxies (incrementally represented by the dashed line) in our 
larger survey of the UDF field. The bottom 4 panels show models of axial ratios for pure- 
clump galaxies with randomly positioned clumps. In each case, there are 5 curves that 
correspond to intrinsic flattenings of Z = 0.1, 0.3, 0.5, 0.7 and 1 (which is an unflattened 
spherical distribution). These curves move progressively to the right in the panels. The 
number of clumps in the models is different for each panel, as indicated. The figure shows 
how UDF galaxies dominated by 5 to 10 clumps, such as chain and clump-cluster galaxies, 
are most likely flattened with an intrinsic axial ratio of 0.1 to 0.3, as shown by the black and 
red curves (in the electronic edition) in the second and third panels up from the bottom. The 
decrease in the distribution at high axial ratio is from the intrinsically irregular structure of 
the galaxies (smooth disks would not have such a decrease), while the decrease at low axial 
ratio is from the minimum intrinsic flattening. Evidently, the clump-cluster galaxies are in 
the process of forming thick disks. 



